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Abstract. Tight Bell inequalities are facets of Pitowsky's correlation polytope and 
are usually obtained from its extreme points by solving the hull problem. Here we 
present an alternative method based on a combination of algebraic results on extensions 
of measures and variable elimination methods, e.g., the Fourier-Motzkin method. Our 
method is shown to overcome some of the computational difficulties associated with 
the hull problem in some non-trivial cases. Moreover, it provides an explanation for 
the arising of only a finite number of families of Bell inequalities in measurement 
scenarios where one experimenter can choose between an arbitrary number of different 
measurements. 

PACS numbers: 03.65.Ud 
1. Introduction 

Bell inequalities [lj play a fundamental role in the investigation of non-classical features 
of quantum mechanics (QM), from the foundational problems raised by the impossibility 
of a local hidden variable description of quantum systems, to applications for quantum 
information processing like entanglement-based reduction of classical communication 
complexity [21 E], device-independent secure communication [H [5] and randomness 
expansion [6]. Similar inequalities, called non-contextuality inequalities, arise in the 
discussion of quantum contextuality [H [8j [9] . 

In fact, the origin of all these inequalities is rooted in logic and probability 
[TP] Hi] 112] and can be traced back [13] to Boole's notion of conditions of possible 
experience |14j ; they are necessary conditions for the interpretation of a set of 
experimentally observed relative frequencies as probabilities in a single probability space. 

From an algebraic point of view [15], Bell and non-contextuality inequalities can 
be seen as conditions for the extension of a function, defined on a subset of a Boolean 
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algebra, to a normalized measure on the entire Boolean algebra, a problem already 
investigated by Horn and Tarski [16]. 

A general method for the derivation of Bell and non-contextuality inequalities is 
based on Pitowsky's notion of correlation polytope [IT] . Such a polytope appears in 
many fields and with different names (see [T7] for a extensive survey of the topic). Its 
geometrical structure is completely defined by a set together with a family of its subsets. 

In the framework of QM, the above set is a set of observables and the subsets 
correspond to subsets of compatible or jointly measurable observables. The operational 
definition of such notions is still under discussion [18J. They have been defined in 
different ways according to the experimental scenario and the physical assumptions 
involved (e.g., space-like separated measurements, sequential measurements). Such 
distinctions are irrelevant for our discussion. On the contrary, the use of an abstract 
notion of compatibility allows us for a unified presentation of our mathematical results 
which are valid in all the above mentioned scenarios. 

Given a set of observables together with the subsets of compatible observables, the 
corresponding correlation polytope is defined as the convex hull of a set of vectors, the 
vertices or extreme points of the polytope, representing possible truth assignments for 
observables and logical conjunctions between compatible pairs, triples, etc. 

Every convex polytope has two representations: One as the convex hull of its 
vertices (V-representation) and the other as the intersection of a finite number of half- 
spaces (^-representation), each given by a linear inequality. Complete sets of tight (i.e., 
facet supporting) Bell inequalities associated with a given set of observables precisely 
amount to the "H-representations of the corresponding correlation polytopes. 

The "H-representation of a convex polytope can be computed starting from its V- 
representation by solving the hull problem. For high-dimensional polytopes, this is a 
very difficult problem. In fact, this method has been used only to compute simple cases 
[HUEO], by means of programs such as cdd [21], Irs [22J and porta [23J. 

More recently, an alternative approach has been proposed by Avis, Imai, Ito and 
Sasaki [2H [25J, which is based on the relation between correlation polytopes and a 
family of convex polytopes called cut polytopes, well studied objects in polyhedral 
combinatorics, and a variable elimination method derived from Fourier-Motzkin method 
of elimination of variables (see, e.g., [26]) and called triangular elimination. Their 
method provided more than two hundred millions new tight Bell inequalities for cases 
in which the corresponding hull problem is computationally intractable. However, since 
such a method is based on cut polytopes, it can be applied only to the bipartite case. 

Here we present an alternative general method based on the application of Fourier- 
Motzkin variable elimination method to conditions derived in Ref. [15j as consistency 
conditions for putting together partial extensions of quantum probabilities in order to 
obtain a classical probability description. More precisely, such conditions are expressed 
in terms of a systems of linear inequalities where also correlations between incompatible 
observables appear as variables: A classical probability space representation exists for a 
given set of QM predictions if and only if the corresponding system of linear inequalities 
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admits a solution; Bell, or non-contextuality, inequalities are obtained by eliminating 
the variables associated to correlations between incompatible observables. Our approach 
can be seen as a generalization of Fine's derivation of the Clauser-Horne-Shimony-Holt 
(CHSH) polytope [10]. 

The above method also provides a generalization of the result obtained by Sliwa (27] 
and Collins and Gisin [20] , namely, the appearance of only a finite number of families of 
Bell inequalities in measurement scenarios where one experimenter is allowed to choose 
between an arbitrary number of different measurements. 

The variable elimination method can also be applied to a family of correlation 
polytopes called complete probability polytopes. In this case, the method may not 
provide a significant computational advantage with respect to other known algorithms. 
Nevertheless, we believe that such results are interesting and we collected them in 
Appendix A 

The paper is organized as follows: In Sec. [2j the notion of correlation polytope 
is reviewed, together with the result on extensions of probabilities associated with tree 
graphs. In Sec. [3j we discuss their application to the problem of the computation of 
^-representation for general correlation polytopes, and provide some examples. Finally, 
in Sec. [| we discuss some possible further developments and computational results. 

2. Correlation polytopes and extensions of measures 

In this section, we recall the definition of correlation polytope and its probabilistic 
interpretation; then we introduce the notion of compatibility graph for a set 
of observables and present the main result about the classical probability space 
representation for the corresponding QM predictions. 

We use the definition of correlation polytope given in [T7J, which is the natural 
generalization of Pitowsky's notion [12] to higher order correlations, but with a slightly 
different notation. 

Given a set of propositions Q = (Ax,...,A n ) and a family X of subsets of Q, 
i.e., X C 2 s , we define the sets Sk, k = 2,...,m with m < n, as the sets of 
logical conjunctions Sk = {A^ A • • ■ A A ik \ ij ^ if, {A^, . . . , A ik } el}. The lines of 
the truth table associated to the above set of propositions and logical conjunctions 
between them, namely the 2 n vectors of KM+l s 2|+-+|S' m | ; 



where e = (£%,..., s n ) G {0,1}", are called the vertices of the correlation polytope. 
Their convex hull, i.e., the set of points generated by their convex combinations, is 
called the correlation polytope associated to X and denoted as COFP{X). 

It is convenient to introduce the following notation which makes apparent the 
correspondence between coordinates and joint probabilities. The coordinates of a point 
p e Rl0l+|Sa|+...+|Sr»| w m be denoted as 




(1) 



V = (Pi, ■■■,Pn,- --Pij, ■ ■ -,Ph...i 



...). 



(2) 
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In order to clarify our notation, it is convenient to present a simple example, 
the CHSH scenario. Consider the bipartite scenario in which one experimenter, say 
Alice, can choose between two dichotomic measurements, associated with propositions 
A\ and A2, and another experimenter, say Bob, can choose between two dichotomic 
measurements, associated with propositions A 3 and A4. The set of propositions is, 
therefore, Q = (Ai, A2, A3, A4). Moreover, the measurements associated with A4 and Aj, 
for i — 1, 2 and j = 3, 4, can be performed jointly and, consequently, it makes sense to 
consider the following set of logical conjunctions S2 = {A1AA3, A1AA4, A 2 AA 3 , A2AA4}. 
The associated polytope is described by 2 4 = 16 vertices in IR 8 , namely 

U £ = (Si, 62, £3, £4, ^1^3; ^1^4; ^2^3? ^2^4); £% G {0, 1}, (3) 

where £i represents a classical {0, l}-valued assignment to proposition Ai and £{£j the 
classical assignment for the logical conjunction Ai A Aj. 

As a consequence of Weyl-Minkowski theorem (see, e.g., p2]), each convex polytope 
has a double description: One as the convex hull of its vertices u E , i.e., the V- 
representation, and one as a (finite) intersection of half-spaces which generates it, each 
one given by a linear inequality, i.e., the ^-representation. 

The convex hull of vertices ^ gives a set of linear inequalities constraining the 
coordinates of a generic point in IR 8 , which we denote by 

V = (Pl,P2,P3,P4,Pl3,Pl4,P23,P24). (4) 



Such inequalities, known as CHSH inequalities, are presented in eqs. (7a)— (7 d) below 
and will be discussed in more detail. The interpretation of such a geometrical object 
is the following: Given a vector belonging to the polytope, each component represents 
the joint probability for the corresponding subset of propositions, e.g., pi$ represents 
the joint probability Prob(Ai A A3); the whole polytope gives possible ranges for such 
joint probabilities if the underlying probabilistic structure is assumed to be classical, 
i.e., given by a probability space or, equivalently (since the number of propositions is 
finite) by a normalized measure on a Boolean algebra. 

It is interesting to define the complete probability polytope associated to 
n propositions Q = {Ai, . . . , A n }, it is the correlation polytope COR n (2 g ), i.e., 
the correlation polytope for a set of compatible proposition, in particular, 

|S| + E*l^l = ELi© = 2"-i. 

The importance of such polytopes has been already recognized in [17], namely, 
the fact that each COR°(l) can be obtained as a projection of COR n {2 5 ) (on the 
subspace IR' 1 '). A more detailed discussion, containing also an algebraic derivation of 



the ^-representation for complete probability polytopes, can be found in Appendix A 



We now recall a result on extension of probability measures which allows for 
another projection-based approach to the computation of ^-representation of correlation 
polytope. More details can be found in Ref . [T5] . 

Proposition 1. (Representation of tree graphs) Consider a set of probabilities Pi on 
a set of {0, 1} -valued observables Ai and correlations pij on a subset of pairs Ai,Aj, 
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defining a probability on each pair Ai,Aj, with pi = (Ai), pj = (Aj), pij = (AiAj); now 
depict observables Ai as vertices and the above pairs as edges in a graph. Then any set 
of predictions associated with a tree graph, i.e., a graph without closed loops, admits a 
classical representation. 

Proposition 2. The same holds with yes/no observables Ai substituted by free Boolean 
algebras Ai, Pi by probabilities on Ai, p^ by probabilities on the Boolean algebra freely 
generated by the union of the sets of generators of Ai and Aj. 

We recall that a Boolean algebra is freely generated by n generators Bi , . . . , B n if 
such generators are as much unconstrained as possible, i.e., they satisfy no conditions 
except those necessary conditions defining a Boolean algebra (e.g., distributive law). 
Since all Boolean algebras freely generated by n < oo generators are isomorphic, 
for the sake of simplicity, we can think of the the algebra of subsets 2 X of 
the set X = {0, l} n , with set theoretic operations (fl,U, c ); then the subsets 
Bi = {(xi, . . . , x n ) G X | Xi = 1 }, for i = 1, . . . , n, can be taken as free generators. 
In terms of propositions and truth assignments, the free Boolean algebra assumption 
amounts to the assumption that each possible {0, l}-valued assignment to propositions 
is admissible. For more details see [2H]- Notice that this is precisely the way how we 
define the 2 n vertices of the correlation polytope. 

In the next section, we shall show how to exploit the above results to obtain 
"H-representation for general correlation polytopes starting from lower dimensional 
polytopes and using Fourier-Motzkin method. 

3. Derivation of Bell inequalities 

The above results provide a general method for the computation of "H-representation for 
correlation polytopes which consists in solving the hull problem for a smaller polytope, 
then constructing an higher dimensional polytope from that solution, and applying 
variable elimination methods such as, e.g., the Fourier-Motzkin method. 

For the convenience of the reader we recall briefly the Fourier-Motzkin method, for 
more details see Consider a system of inequalities of the form Ax < b, where A 
is a m x d real matrix, x = (x\, . . . ,Xd) £ K d and b 6 IR m , and suppose we want to 
eliminate the variable Xd- After a proper normalization, each inequality will be of the 
form OjiXi-K . .-t-a^d < h with a i( i being +1, —1 or 0. By summing an inequality where 
did = +1 with an inequality where ajd = —1, we obtain a new inequality not containing 
Xd- By repeating the above operation for every possible pair i, j such that = +1 and 
ajd = —1, and considering also the inequalities where = 0, we obtain a new system 
of inequalities not containing Xd- From a geometric point of view, since the system of 
linear inequalities can represent a polytope (more generally, a cone), the above operation 
amounts to a projection on the coordinates associated with the variables x±, . . . , Xd-i- 

Our method consists in exploiting the automatically existing classical representa- 
tions for subsets of observables with compatibility relations described by tree graphs, see 
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(a) (b) (c) {d) 

Figure 1. (a) Graph of the compatibility relations between the observables in 
the CHSH scenario, (b) Partition in two subset of three observables with intersection 
{A\, A2}. (c) Tree graph obtained by extending the probability measure on the algebra 
generated by Ax^A^. (d) Asymmetric case with additional observables on Bob's side. 

Propositions [l] and [2] above. Conditions for classical representability arise as consistency 
(i.e., coincidence on intersections) conditions for putting together partial extensions as- 
sociated with subgraphs, giving rise to a description of the initial compatibility graph as 
a tree graph on such extended nodes. Such consistency conditions are expressed in terms 
of the existence of a solution for a set of linear inequalities. One of the main application 
of Fourier-Motzkin algorithm is precisely deciding whether a system of inequalities has 
a solution. 

Moreover, we shall discuss how this result gives an account for the appearance of 
a finite number of families of Bell inequalities in cases in which one experimenter can 
choose between arbitrary number of different measurements while such a number is 
fixed for the others, as recognized by Sliwa [27] and Collins and Gisin (20] for the case 
in which Alice perform two measurement and Bob an arbitrary number. 

We recall that, although our method is general, different strategies are possible 
corresponding to different partitions of the initial graph into subgraphs. We shall discuss 
our method by means of some simple examples. The first one is the derivation of the 
CHSH polytope. It is interesting to notice that it is analogous to that presented by Fine 
[TO] ; our method can be seen as a generalization of his idea to an arbitrary number of 
observables. 



3.1. CHSH polytope from Bell-Wigner polytope 

The CHSH polytope is generated by the following set of vertices 

u £ = (ei, 52, £3, £4, £i£3, £i£4, £2^3, £2^4), £ {0,1}. (5) 

It is associated with a bipartite measurement scenario in which Alice can choose between 
two measurements, associated with propositions A\ and A 2 , and Bob can choose between 
two measurements, associated with propositions A 3 and A4. 
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As already recognized by Fine \±0\ and discussed also in [151 EE], the existence 
of a classical description for the four observables is equivalent to the existence of 
classical descriptions for the two subsystems, {Ai,A 2 , A 3 } and {A±, A 2 , A4}, coinciding 
on {A±, A2}. In fact, the two classical descriptions would give rise to an extension of the 
probability assignment to the four observables satisfying the hypothesis of Proposition 
I (see Figure [l] (a),(b),(c)). 

The constraints on the subsystem {Ai, A 2 } imposed by the third observable, A 3 or 
A4, are described by the Bell-Wigner polytope, i.e., the correlation polytope associated 
to three propositions and their pairwise logical conjunctions, which is given by the 



following inequalities, obtained in Ref . [T2] . 

< Pij <pi, < pij < pj, ij = 12, Is, 2s, (6a) 
Pi + Pj ~ Pij < 1, ij = 12, Is, 2s, (66) 

1 ~Pl -P2 ~Ps +P12 +Pls +P2s > 0, (6c) 

Pi - P12 - Pis + P2s > 0, (6d) 

P2 ~ P12 ~ P2s + Pu > 0, (6e) 

Ps -Pis ~P2s +P12 > 0, (6f) 



for s = 3 or 4. For fixed s (3 or 4), the above inequalities give ranges for the classical 
probabilities for a system of three proposition and their pairwise joint probabilities. 
It is interesting to notice that, by Proposition [TJ each subset of three observables 
always admit a classical representation (see also Figure [I] (b)); therefore, for fixed s, 
the above system of inequalities admits a solution, i.e., a value can be assigned to the 
joint probability p\ 2 (which is assumed to be experimentally unmeasurable, since it 
is associated with two noncommuting observables) consistently with the constraints 
imposed by the classical model and the measurable joint probabilities p\ s and p 2s . 
This implies that, for fixed s, the elimination of variable p\ 2 gives rise only to trivial 
inequalities, i.e., inequalities always satisfied by quantum predictions. An analysis of 
the constraints imposed on such "experimentally unmeasurable" joint probabilities has 
been presented in |29j. Analogous situations always arise in the bipartite scenario 
and this additional property just gives more information about the triviality of certain 
inequalities. However our method can also be applied to the general multipartite case. 

Classical representability for QM prediction in the CHSH scenario amounts, 
therefore, to the existence of a common solution for the two systems of linear inequalities, 
namely, it amounts to the existence of a value for p\ 2 consistent with the constraints 
imposed by classical descriptions for the two subsystems of three observables. 

As a consequence of general properties of Fourier-Motzkin method (see [26]), a 
system of inequalities admits a solution if and only if the projected system, i.e., the 
system obtained by eliminating one or more variables, admits a solution. It follows that 
the above set of measurements admits a classical description if and only if measured 
correlations, i.e., correlations between compatible observables, satisfy the system of 
inequalities obtained by eliminating p 12 . 
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In order to eliminate P12, we just combine inequalities where P12 appears with 
opposite sign and keep inequalities where it does not appear. As discussed before, by 
Proposition [TJ combining inequalities with the same index s gives rise only to trivial 
inequalities, in particular, they are also redundant. A posteriori, we know that the same 
happens also for inequalities (6a) and (66). 



The interesting inequalities are those obtained from the combination of (6c)-(6j) 
for different s, namely 



- 1 < P13 + Pu + P24 ~ P23 - Pi ~ Pa < 0, 

-1 < P23 + P24 + Pl4 ~ Pl3 ~ P2 ~ Pi < 0, 
-1 < PU + Pl3 + P23 - P2A ~ Pi ~ P3 < 0, 
-1 < P24 + P23 + Pl3 - PU ~P2-P3<0. 



(7a) 
(7b) 
(7c) 
(7d) 



For instance, — 1 < p\ 3 + pu +P24 —P23 —Pi~Pa is obtained as a sum of inequality (6c) 
for s = 4, which contains the term +P12, with inequality (6e) for s = 3, which contains 
the term — p 12 . 



Inequalities (7a)-(7d), together with the inequalities ( |6a[ ), (6b) in which p i2 does 
not appear, give the "H- representation of the CHSH polytope (compare with, e.g., [T2]). 



3.2. Bipartite (2,n) scenario 

An analogous argument applies to the scenario in which Alice can choose between 
two measurements and Bob can choose among n > 2 measurements, associated to 



propositions A 3 , . . . , A n+2 : The initial system of inequalities is still given by (6a)-(6/), 
but with s taking values in {3, 4, . . . , n + 2} (see Figure [l] (d)). 

Since only one variable (i.e., P12) has to be eliminated, at most two inequalities with 
different index s can be combined to give a valid inequality. Therefore, the final set of 



inequalities is given by (7a)-(7d) with the pair 3,4 substituted by any pair i,j with 



i, j G {3, . . . ,n + 2} and % < j. This is precisely the result obtained in Refs. [20l [27]. 



3.3. Two parties, three settings 

Now consider the bipartite scenario in which Alice can choose among three 
measurements, associated with propositions A\,A2 and A3, and Bob can choose among 
three measurements, associated with propositions A 4 , A 5 and A 6 . 

Analogously to the previous case, see Figure |2j the existence of a classical 
description for the six observables is equivalent to the existence of classical descriptions 
for the three subsystems, {At, A 2 , A 3 , A^}, {A\, A 2l A3, A$} and {Ax, A 2 , A 3 , A 6 }, 
coinciding on {Ai, A 2) A 3 }. A probability on {Ai, A 2 , A 3 } is completely defined, see 
Lemma [5] below, once the probabilities Pi, p 2 , P3, P12, P13, P23, P123 are given. 

It is therefore sufficient to calculate the correlation polytope associated with 
probabilities Pi,P2,P3,Ps,Pis,P2s,P3s,Pi2,Pi3,P23,Pi23, then consider the system given 
by all the above inequalities for s = 4, 5, 6 and eliminate the variables pi2,Pi3,P23,Pi23- 
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Figure 2. (a) Graph of the compatibility relations between the observables in the 
bipartite (3,3) scenario, (b) Tree graph obtained by extending the measure on the 
algebra generated by A\, A 2 , A3. 

3.4- Bipartite (3, n) scenario 

Again, by adding observables only on Bob's side, one just obtains more copies of the 
initial system of inequalities, but with different indices s. The situation is analogous to 
that depicted in Figure [I] (d), but with A12 substituted by A123 

In particular, since four variables have to be eliminated, at most 2 4 = 16 inequalities 
with different indices s can be combined. As a result, all families of valid inequalities for 
the general case (3,n) already arise in the case in which Alice performs 3 measurement 
and Bob 16. 

3.5. Multipartite (m, . . . , m, n) scenario 

The above argument can be extended to the case of p parties in which the first p — 1 
can choose among m, with m fixed, measurements, while the last one can choose among 
n > m, with n arbitrary, measurements: All families of inequalities can be obtained by 
studying the case in which the last experimenter performs 2 k measurements, where k is 
the number of variables to be eliminated. 

4. Conclusions and computational results 

We have presented an alternative method for the computation of half-space 
representation for correlation polytopes based on algebraic conditions and variable 
elimination. A reasonable question is: Does it provide any advantage with respect 
to existing methods? In order to show the advantages of the tree graph method, 
we have computed the ^-representation for some simple polytopes, with (i) our tree 
graph method using existing software implementing the Fourier-Motzkin algorithm; 
specifically, we used porta [23] and FM library [30], and (ii) using standard software 
for solving the hull problem; specifically, we used cdd. 

For simple cases like the (2, 2) (i.e., the CHSH) and (3, 3) scenarios, the computation 
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is equally fast with both methods. However, remarkably, our tree graph method is 
noticeably faster to compute asymmetric scenarios: 

For the (3,4) scenario, the tree graph method implemented with porta completed 
the calculation in rj 11 minutes, while cdd needed ~ 20 minutes. The 11 minutes include 



the time (seconds) required to calculate the initial polytope (see Sec. 3.3). 

For the (3, 5) scenario, the tree graph method implemented with porta completed 
the calculation in r* 72 minutes, while cdd was still running after a week and we had 
to stop it. All computations were performed on the same machine with an Intel Xeon 
CPU running at 3.20 GHz. 

However, beyond the practical advantage the tree graph method has to compute 
some scenarios, the main purpose of this paper is to point out that, for some polytopes, 
the ^-representation can be easily expressed by means of simple algebraic arguments, 
and then these polytopes can be related to those associated to polytopes arising in 
physically interesting scenarios by means of variable elimination methods. For example, 
the Boolean algebraic approach in Ref. [E] provided, already in 1948, complete sets of 
Bell inequalities. To our knowledge, this approach has never been investigated in relation 
to Bell inequalities, and we hope that our work will stimulate further developments in 
this direction. 
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Appendix A. 

In this appendix we recall some results, Proposition [4] and Lemmas [3j |5j arising from 
the translation of Pitowsky's notion of correlation polytope into the Boolean framework, 
and use them to derive the ^-representation for the complete correlation polytope. For 
more details on such results and the basic notions of Boolean algebra involved, see [TS] . 

After completing our work, we realized that i) obtaining a correlation polytope as 
a projection of the corresponding complete probability polytope need approximatively 
the same amount of time and resources as the usual convex hull algorithms based on 
Fourier-Motzkin elimination, and, ii) similar results were already known in the field of 
integer 0/1 programming. In fact, in Ref. |31j. the ^-representation of the complete 
probability polytope is obtained by finding the explicit form of the matrix N we defined 
in Lemma [9] below. Nevertheless, we believe our approach is interesting, since it is 
based on logical and probabilistic, rather than geometrical, notions and provides an 
alternative point of view which has not been explored, to our knowledge, in the context 
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of Bell inequalities. 

First we present the "H-representation for the complete probability polytope 
associated with n propositions. It is given by the following set of 2 n tight inequalities: 

h(ei,...,e n )>0, £ = (£!,..., s n ) e {0,1}™, (A.l) 

where h(si, . . . , £ n ) is defined as 

s k 

h(l, . . . , 1,0, . . . ,0) = PL..S - ^2 Pi-si + Pi-difaH 

s + l<l<s + k s + l<; 1 <i 2 <s + fe (J± 2) 

•••+(-i) fe - 1 Yl pi-shi 2 -i^ + (-i)Vi..., 



■in 

hi<;i<(2<— <ik-i<s+k 



for all 1 < s < n and s + A; = n; all other values, except /z(0, . . . , 0), are defined by a 
permutation of indices {1, . . . , n} providing the above form, i.e., (1, . . . , 1, 0, . . . , 0), for 
(ei, . . . , £ n ); in the case s = 0, the first term on the r.h.s is 1. 

The notion of correlation polytope has been introduced by Pitowsky [12J in terms of 
propositional logic, however the same problem can be expressed in terms of free Boolean 
algebras with the identification (see of atomic proposition with free generators, logic 
operations with Boolean operations, truth assignments with two-valued measures and 
probability assignments with normalized measures. 

To prove that the above representation is correct, we first need the following basic 

fact 

Lemma 3. Let *B be a finite Boolean algebra with k atoms {a 1; . . . , a^}; then a function 
f : X C 53 — > [0, 1] can be extended to normalized measure fi on 03 if and only if there 
are k numbers Ai, . . . , A^, with Aj > and Yli=o ^ = ^ such that 

k 
i=l 

where 5 ai is the two-valued measure which is 1 on a^. /i is given by 
k 

// = ^V a< , (A.4) 

i=l 

and the coefficients correspond to measure of atoms, i.e., Aj = /i(aj). 

As a consequence, Pitowsky's result can be expressed as follows: 

Proposition 4. Let 03 be a Boolean algebra freely generated by Q = {Ai,...,A n } 
and consider X C 23 with X = {Ai, . . . , A n , . . . , A^ D Aj, . . . , A^ D Aj H Ak, ■ ■ ■}, i.e., 
X = Q U S2 U . . . S m , m < n, where elements of Si are the intersections of I distinct 
generators, but not necessarily all of those possible, i.e., \Si\ < m. Now consider 
f : X — > [0, 1] and define the vector 

P = (pi, • • • ,Pn, ■ --Pij, ■ ■ -,Ph...i m , ■ ■ •) e (A.5) 
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which has as components the values assumed by f on X , namely 

Pi =f(Ai), i = l,...,n, (A.6) 

Pij = f{A i nA j ), Ai PI Aj G 5*2, . . . (A.7) 

Pl ,., ra = /(4n...nA i j, 4n...n4^ m . (A.8) 

For every e G {0, l} n define the vector u e G {0, given 6?/ 

M e (^1; ■ ■ ■ j ^n; ■ ■ ■ j E-i^ji ■ ■ ■ i ^ii^i2 ■ ■ ■ ^i m i •••)■> (A-9) 



i.e., for every component p^..,^ of p there is a corresponding component of u e given by 

£ h ■ ■ ■ ^ik ■ 

Then f can be extended to a normalized measure on 23 if and only if there are 2 n 
numbers A(e), e G {0, l} n , such that 

e£{0,l} n 

with 

A(e) > Ve 6 {0, l} n , and X ( £ ) = 1 ■ ( A - n ) 

ee{0,l} n 

The convex hull of the vectors u £ defined above is known as the correlation polytope; 
clearly, each Ai is associated to a QM observable, each set Sk to a joint measurements of 
k observables, and numbers Pi lv ..,j fc to probabilities for the corresponding measurements 

IfX contains all possible intersections between generators, i.e., \X\ = Ylk=i u) = 2 n — 1, 
it will be called a complete set of measurements and the corresponding polytope will 
be called complete probability polytope, and it will be denoted as CV n ; such names 
come from the fact that there is a bijection between the points of the polytope and 
the normalized (i.e., probability) measures on the Boolean algebra freely generated by 
Ai, . . . ,A n , as a consequence of the following general fact [15]: 

Lemma 5. Let 03 be a Boolean algebra freely generated by {Ax, . . . , A n }, let \x be a 
normalized measure on 05 and let X C 23 be a complete set of measurements. Then the 
measure /i is uniquely defined by the values it assumes on the set X . 

Now consider a Boolean algebra 53 freely generated by {Ax, . . . ,A n }, a complete 
set of measurement X C 53 and a function / : X — > [0, 1]. If / admits extension 
to a normalized measure fi, then the measure of the atoms of 23, i.e., the elements 
A\ l n A £ 2 2 n . . . n A £ n n , where G {0, 1} and Af is defined as A t if = 1 and A\ (the 
Boolean complement) otherwise, can be written as 

n(Ax n . . . n AO = f(A x n . . . n A n ) , 

n{Ax n . . . n A c n ) = f{Ax n . . . n A„_0 - f{A x n...ni„), 

! : ! (A.12) 

n{A\ n . . . n AD = l - /(A) + Ei<j /(4 n A,) 

- J2 i<j<k /(A n A n A) + . . . + (-i)"/(A n...n A») , 
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where the above expressions are obtained by using the identities fi(A) = fi(A D B) + 
fi(A PI B c ) for all A, B e 03 (additivity of measure on disjoint elements) and /i(l) = 1, 
1 being the identity in 03 (normalization). 

The exact expression for the values of the measure on atoms in terms of values on 
intersections is given by the following general result: 

Lemma 6. Given a Boolean algebra 03 and n > 2 elements Bi, . . . , B n e 03 and a 
measure /x on 03, it holds 

ii{B x n...nfi s n 5 S C +1 n . . . n B c s+k ) = ^{B, n . . . n £ s )+ 

- /i(5 1 n...nB s nB i )+ ^ ^(s 1 n...ns a ns, 1 ns, 3 ) + ... 

s + l<Z<s + fc s + l<Z 1 <Z 2 <s + fc 

( A. lo J 

■•• + (-1) ^ n . . . n s s n B h n B la n . . . n B lk _J+ 

s+i<i 1 <; 2 <...<i fc _ 1 <s+fe 

+(-i)V(£in...nfl a+fc ) , 

TOi/i s and sttc/i t/iat 1 < s < n, and fc = n — s. In the case s = 0, i.e., ^(B^D. . ,C\B%), 
an analogous expression holds, but with the first term fi(Bi Pi ... Pi B s ) substituted by 
/i(l), which is equal to 1 if fi is a normalized measure. 

Proof Fixed n > 2, we proceed by induction on k. For k = we have nothing to 
show. For k — 1 it is sufficient to use the identity fi(A fl B c ) = fi(A) — fi(A fl B). 

For the inductive step, it is sufficient to define B[ = B\ fl B c s+k1 apply the inductive 
hypothesis to fJ,(B[ fl B 2 fl . . . fl B s fl B c s+1 fl . . . fl Bg +k _^) and then apply again the 
identity fx(A n B c ) = fi(A) - fj,(A H 5) to eliminate B c s+k . 

It is obvious that the above argument does not depend on n. The proof for s = 
is analogous. ■ 

Now consider a generic function / : X — > [0, 1] and define X as the set of atoms 
X = {Af n# 2 2 n...nA^» | s { e {0, 1} } and a function / : X — > E in terms of / 
in analogy with the result of Lemma [HJ i.e., for each A^ 1 fl fl . . . fl A 1 ^ 1 , f is defined 



as the r.h.s. of the corresponding equation in (A. 12). Notice that by construction and 
independently of the values of /, in particular, independently of the existence of an 
extension for /, / satisfies 

f(A h n . . . n A ik ) = Yl Mi n • • • n A, n A^ n...n^), (A.14) 

£i k+l >—' £ in 

1 = J2 f(A?n...nA?) , (A.15) 

El,...,£ n 

where (ii, . . . , i n ) is a permutation of (1, . . . , n). In particular, each of the two functions, 
/ and /, is completely determined by the other one. 
Extensions of / and / are related by the following 

Lemma 7. / admits an extension to a measure •<=>■ / admits an extension to a measure. 
Moreover, if the two extensions exist, they coincide. 
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Proof The first implication follows from the definition of /; the other follows from 



Eq. (A. 14) together with the identity 

A h n...nA ik = |J ^n...n4n^Vn...n^, (A.16) 

for all (ix, . . . , i n ) permutation of (1, . . . ,n) and all 1 < k < n, since the above is a 
disjoint union. 

To conclude, it is sufficient to notice that a measure is completely defined by the 
values it assumes on the atoms of the Boolean algebra (since each element can be written 
in a unique way as as disjoint union of atoms) and that each of the two functions is 
uniquely determined by the other one. ■ 

We are now able to prove the main result 

Theorem 8. (H -representation of CV n ) 



A vector p = (p h . . . ,p n , . . . ,p y , . . . ,p h ... im , . . . ,pi. 
components satisfy the following 2 n inequalities 

< pi,.. n 

< Pl...( n -1) — Pl...(n-l)n 



< 1 - Ei Pi + Eii Pij ~ E«fc Vijk + ■■■ 



»2 n -l 



belongs to CV r 



its 



(A.17) 



where the r.h.s. of inequalities (A. 11) are obtained from the r.h.s of eqs. (A. 12) by 
substituting f(A i:L D . . . D A ik ) with p^...^ for all . . . , i k } C {1, . . . , n}. 

Proof Given p e IR 2 ™ -1 , consider the function / : X C 53 — > R, where 25 
is the Boolean algebra freely generated by {Ai,...,A n } and X a complete set of 
measurements, defined as 

f(A H n...f]A ik ) =p h ... ik , for all {i u ...,i k } C {l,...,n} . (A.18) 

Let us assume that / takes values in the interval [0, 1], we shall prove it later. Then, by 
Proposition [4], the above problem amounts to the problem of the existence of a measure 
extending / on 33, and therefore, by Lemma [7j to the problem of the existence of a 
measure extending /. Since / is defined on the atoms of 2$ and each element can be 
written uniquely as a disjoint union of atoms, / can be defined on each B e 55 as the 
sum of the valued assumed on the atoms A^ 1 Pi ... Pi A^ 1 C B. In this way we obtain 
a unique extension of / on 03 which is additive on disjoint elements and normalized by 
construction (see Eq. (A. 15)). It follows that such an extension is a measure <^=^ it 
is non-negative <^=^ / is non-negative on the elements of X. By construction of /, see 



Eq. (A. 12), such a condition is equivalent to the system of inequalities (A.17). 



To conclude, we just have to show that the function / defined above takes values 
in [0,1]. Since inequalities (A.17) amount to the condition / > 0, by Eq. (A. 14), it 



by Eq. (A.15), also that / < 1. 



follows that / > 0. This implies, 

We now prove that the above inequalities are tight, i.e., they cannot be represented 
as a positive weighted sum of other inequalities and define facets of the polytope. First, 
we need the following 
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Lemma 9. CV n is full dimensional in M. 2 " 1 for all n. 

Proof Given p G CV n , it can be written, with the notation of Proposition |4j 
as P — J2 £ ^( £ ) u e- Now, consider the vectors p',u' £ ,X G IR 2 " defined as p' = (p, 1), 
u' e = (u £ , 1) and A as the vector given by the 2 n values A(e) for e G {0, 1}™. Notice that, 
by Lemma p[ the vector A is given by the values assumed by / 



Given p G CV n , by definition of CV n and Eq. (A. 15), there exists a 2™ x 2™ matrix 
M, which columns are given by the vectors u' £ , such that p' = MX. On the other 
hand, by definition of /, each component of A can be written as a linear combination 
of the components of p', i.e., there exists a matrix iV such that A = Np'. Since the set 
of admissible values for A, see proof of Theorem [sj contains a basis for IR 2 ™, e.g., the 
canonical basis, it follows that iV = M~ l . 

As a consequence, the columns of M, i.e., the vectors u' £ , are linearly independent. 
This implies that the dimension of the subspace spanned by u £ is 2 n — 1. ■ 

In order to show that the inequalities (A. 17) are tight, it is sufficient to notice that 
CV n has exactly 2 n facets. It is sufficient to count the number of affine hyperplanes, 
i.e., affine subspaces of dimension 2™ — 2, generated by the vertices u £ . By the above 
lemma, the vectors u £ are affinely independent, therefore each subset of 2 n — 1 vectors 
defines an affine hyperplane. Therefore, the number of facets of CV n is ( 2 ™_i) = 2 n . 
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